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Abstract 

We discuss certain aspects of the production of ®*^Fe in massive stars in the range 
between 11 and 120 Mq, both in the hydrostatic and explosive stages. We also 
compare the ®°Fe/^®Al 7 -ray line flux ratio obtained in the present calculations to 
the detected value reported by INTEGRAL/SPI. 
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1 Introduction 


The nucleus ®°Fe is a long lived (r ~ 2 Myr) radioactive isotope that should 
be present in an appreciable amount in our galaxy. Historically, ®°Fe has been 
considered as a key isotope to understand whether or not massive stars are 
the main contributors to the diffuse ^®A1 present in the Galaxy, another long 
lived radioactive isotope (r ~ 1 Myr) traced by the detected 1.809 MeV y-ray 


emission line at a 
80’s 


evel of ~ 4 X 10 ^ cm ^ s Indeed, already in the early 


ClavtonI ()1982I1 pointed out that SNII are the only candidate sources for 


^®A1 to produce also ®°Fe and hence that the detection of this isotope in the 
Milky Way could constitute a strong argument in favour of SNII as the main 
contributors to the galactic ^®A1. The first detection of ^ °Fe in the Galaxy 
was obtained with RHESSI and reported bv ISmithI (j2003ll . The line flux de¬ 
tected implies a ^°Fe/^^Al y-rav l ine flux ratio of 0.16 (for each ®°Fe line). 


More recently 


Harris et ah 


( 200511 reported the hrst detection of °°Fe decay 


lines at 1.173 MeV and 1.333 MeV in spectra taken by the SPI spectrom¬ 
eter on board INTEGRAL during its first year, yielding a y-ray line flux of 
3.7±1.1 X 10“® y cm“^ s“^. The same analysis applied to the 1.809 MeV line of 
^®A1 yielded a ®°Fe/^®Al y-ray line flux ratio of ~ 0.11T0.03. From a theoretical 
side, many groups have performed calculations of nucleosynthesis in massive 
stars, estimating the amounts of either ^®A1 or both ^®A1 and ®°Fe ejected in 
the interstellar medium. However, no set of models covers an extended grid of 
stellar masses. Indeed, several research groups interested in the presupernova 
evolution and explosion of ms 
for stars up to a mass of 40 


2mm 


Wooslev fc Weaver. 






stars provided vields of be 

)th ^®A1 and ®°Fe 

Ghieffi & Limoneri. 

2W4; 

Rauscher et ah. 


Thielemann et al 


199611 - among these works 
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only 


RauscheretaL 


( 200211 included mass loss in the computations. On the 


other hand, groups mainly interested in the evolution of massive stars includ¬ 


ing mass loss computed the evolution up to the end o 


Palacios et ah. 

2W5; 

Lanerer et ah. 

1995) 


central He burning 


Mevnet et al 


1997 


In this paper we will discuss to some extent the production of ®°Fe in mas¬ 
sive stars in the range between 11 and 120 Mq, both in the hydrostatic and 
explosive stages. We will also provide theoretical predictions for the ®°Fe/^®Al 
7 -ray line flux ratio of such a generation of massive stars and we will compare 
them with the observations. 


2 The stellar models 

The yields of ^®A1 and ®°Fe discussed in this paper are based on a new set of 
presupernova models and explosions of solar metallicity stars, with mass loss, 
in the mass range between 11 and 120 Mq, covering therefore the full range of 
masses that are expected to give rise to Type Il/Ib/Ic supernovae as well as 
those contributing to the Wolf-Rayet populations. All these models, computed 
by means of the latest version of the FRANEC code, will be presented and 
discussed in a forthcoming paper (Limongi & Chiefh in preparation). 


3 The production of ®°Fe in massive stars 

The isotope ®°Fe is an unstable nucleus (its terrestrial decay time is r ~ 
2 X 10®y) that lies slightly out of the stability valley, its closest stable neigh- 
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Fig. 1. Total yield of as a function of the initial mass (filled stars); the hydro¬ 
static shell He burning contribution (filled circles); the hydrostatic shell C burning 
contribution (filled triangles); the explosive contribution (filled squares). 

bor being ^®Fe. At temperatures below 2-10® K ®°Fe is mainly produced by 
neutron capture on the unstable nucleus ®®Fe and destroyed by the (n, 7 ) re¬ 
action (that always overcomes the beta decay). Since ^®Fe is unstable, the 
®°Fe production rate depends on the competition between the ^®Fe(n, 7 )®°Fe 
reaction and the ^®Fe(/5“)®®Co decay. At temperatures above 2 • 10® K, ®®Fe 
is totally destroyed mainly by ( 7 , n) photodisintegrations and (p, n) reactions. 
Hence, two requirements must be fulhlled in order to obtain a substantial 
production of ®®Fe: 1 ) the neutron density must be high enough (i.e. the tem¬ 
perature must be high enough) to allow the ®®Fe(n, 7 )®®Fe reaction to over¬ 
come the ^®Fe(/5“)®®Co decay; 2 ) the temperature must be not too high to 
lead to the complete destruction of ®®Fe via ( 7 , n) and (p, n) reactions. An 
order of magnitude estimate of the neutron density necessary to cross the 
^®Fe bottleneck may be derived by equating the (n, 7 ) and f3~ decay rates. A 
comparison between such a neutron density and the central neutron density 
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obtained during the evolution shows that no substantial production of ®°Fe is 
obtained in any central burning of a massive star. Indeed, for temperatures 
below 2 ■ 10® K the actual neutron density is always below the one required 
to produce a substantial amount of ®®Fe. On the contrary, the larger burning 
temperatures at which the shell burnings occur allow a much higher produc¬ 
tion of ®®Fe. In particular, in stars in the mass interval 40-120 Mq, shell He 
burning occurs at temperatures as high as 4 • 10® K. This implies a neutron 
density of 6 x 10^® — 10^^ n/cm®, the neutrons being produced mainly by the 
^^Ne(Q;, n)^^Mg reaction, and hence a large amount of ®®Fe. In analogy with 
shell He burning, also shell C burning occurs at temperatures high enough 
(T > 1.3 ■ 10® K for stars in the mass interval 20-120 Mq) that a high neutron 
density is obtained (6 x 10^^ —6 x 10^^ n/cm®) and hence a large amount of ®®Fe 
is synthesized. Also in this case the main neutron source is the ^^Ne(Q;, n)^®Mg 
reaction, the alpha particles being provided by the Q;)^®Ne reaction. 

It must be noted, at this point, that the presence of a convective shell plays a 
crucial role for the synthesis of ®®Fe. Indeed, it has the double responsibility of 
bringing new fuel {a particles and ^^Ne) in the region where the active burning 
occurs and simultaneously of bringing the freshly made ®®Fe outward in mass, 
i.e. at lower temperatures where the neutron density becomes negligible and 
the ®®Fe half life increases substantially. Shell Ne burning occurs at tempera¬ 
tures high enough, but still below the critical value for the total destruction 
of ®®Fe, to allow a large neutron density. However, the lack of an extended 
and stable convective shell lasting up to the moment of the explosion prevents 
the build up of a signihcant amount of ®®Fe. In shell O and Si burnings the 
temperature is so high (above 2 • 10® K) that no appreciable amount of ®®Fe 
is produced. The hydrostatic production of ®®Fe as a function of the stellar 
mass is shown in Figure 1. A comparison between the contribution due to 
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the shell He burning and shell C burning shows that in stars below 60 Mq 
the hydrostatic production of ^‘^Fe is dominated by the contribution of the 
shell C burning while above this limit the ®°Fe is mainly due to the shell He 
burning. The local minimum corresponding to the 30 Mq is the consequence 
of the lack of an efficient C convective shell lasting up to the presupernova 
stage in this model. The isotope ®°Fe is eventually produced during the ex¬ 
plosion in the regions heated up to a temperature of Tpeak ^ 2.2 X 10^ K. 
Indeed, for the typical explosive burning timescales (~ 1 s), below this critical 
temperature the neutron density is too low to allow a substantial production 
of ®°Fe while above this limit ®°Fe is totally destroyed by photodisintegration 
reactions. Since in most cases such a temperature is reached either at the base 


or within the C convective shell and since the matter behind t 


can be assumed radiation dominated (Weaver fc Wooslev 


re shock wave 


198011 . the amount 


of ®°Fe produced during the explosion will depend on the local abundances of 
^°Ne, ^^Na, and ^^Ne left by the last C convective s hell episode and on 


the h nal MR relation at the moment of the core collapse (jChieffi fc Limongi . 
200211 . Figure 1 shows the total yield of ®°Fe as a function of the stellar mass 


together with the explosive and hydrostatic contributions. There is a global 
direct scaling with the initial mass and a quite monotonic behavior: the two 
exceptions are the 13 Mq and the 30 Mq models. The 13 Mq model consti¬ 
tutes a striking exception because a large amount of ^‘^Fe is synthesized by the 
explosion in this case. The reason is that the peak temperature of 2.2 x 10® K 
occurs beyond the outer border of the C convective shell where the abundances 
of and ^^Ne, in particular, are much higher than in the C convective shell. 
The minimum corresponding to the 30 Mq model has been already discussed 
(see above). Below 60 Mq the total yield is dominated by the contribution of 
the C convective shell while above this mass it is the He convective shell to 
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play the major role. The explosive burning almost always plays a minor role. 



Fig. 2. ®^Fe/^®Al 7 -ray line flux ratio integrated over a single power low IMF as a 
function of the slope of the IMF for three different upper mass limits; Mtop = 60 Mq 
(filled triangles); M^p = 80 Mq (filled squares); M^p = 120 Mq (filled circles). The 
solid lines refer to the standard models, while the dashed lines to models in which 
the contribution to the ®*^Fe of the He convective shell of the more massive stars is 
removed. 


4 Comparison to observations 

To compare properly with the observations (see introduction), the yields of 
both ^®A1 and ®°Fe must be integrated over a stellar Initial Mass Function 
(IMF) 0(m). We assume here that the IMF is described by a single power 
low, i.e., (j){m) = Figure 2 shows the ®'^Fe/^®Al 7 -ray line flux ratio 

as a function of the slope x of the IMF, obtained for three different val¬ 
ues of the IMF upper mass limit, i.e., M^p = 60, 80 and 120 Mq (solid 
lines). The horizontal shaded area refers to the ®°Fe/^®Al 7 -ray line flux ra- 
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tio reported by INTEGRAL/SPI ( Harris et al 


2nn5li while the hatched areas 


correspond to the region where the ratio between Type Ib/c and Type II sn- 
pernovae is compatible with the observed valne of 0.3 ± 0.04, as reported by 


fjCappellaro fc Tnratto . 


200lh . Figure 2 clearly shows that our theoretical pre¬ 


dictions always overestimate the observed ®°Fe/^®Al y-ray line flux ratio for 
any choice of the slope of the IMF and the IMF upper mass limit. This is due 
to the copious ®°Fe production in stars more massive than 40 Mq. Indeed, the 
£t is much more improved if such a contribution is artificially removed in these 
stars (dashed lines in Figure 2). A more detailed discussion of the implications 
of these result will be addressed soon in a forthcoming paper. 
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